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^ ■ ABSTRACT 

' We present the results of a deep HST/WFPC2 imaging study of 17 quasars at 2 ~ 

CO , 0.4, designed to determine the properties of their host galaxies. The sample consists 

of quasars with absolute magnitudes in the range — 24 > My > —28, allowing us 
to investigate host galaxy properties across a decade in quasar luminosity, but at a 
single redshift. Our previous imaging studies of AGN hosts have focussed primarily on 
quasars of moderate luminosity, but the most powerful objects in the current sample 
] have powers comparable to the most luminous quasars found at high redshifts. 

QQ ■ We find that the host galaxies of all the radio-loud quasars, and all the radio- 

' quiet quasars in our sample with nuclear luminosities My < —24, are massive bulge- 

CO , dominated galaxies, confirming and extending the trends deduced from our previous 

' studies. From the best-fitting model host galaxies we have estimated spheroid and 

I hence black-hole masses, and the efficiency (with respect to the Eddington luminosity) 

Q^. with which each quasar is emitting radiation. The largest inferred black-hole mass in 

I ' our sample is Mbh — 3 x IO^Mq, comparable to the mass of the black- holes at the 

O I centres of M87 and Cygnus A. We find no evidence for super-Eddington accretion 

rates in even the most luminous objects (0.05 < L/LEdd < 1-0). 
^ ' We investigate the role of scatter in the black-hole: spheroid mass relation in deter- 

mining the ratio of quasar to host-galaxy luminosity, by generating simulated popula- 
tions of quasars lying in hosts with a Schechter mass function. Within the subsample 
of the highest-luminosity quasars, the observed variation in nuclear- host luminosity 
ratio is consistent with being the result of the scatter in the black-hole: spheroid re- 
' lation. Quasars with high nuclear-host luminosity ratios can be explained in terms 

of sub-Eddington accretion rates onto black-holes in the high-mass tail of the black- 
hole:spheroid relation. Our results imply that, owing to the Schechter function cutoff, 
host mass should not continue to increase linearly with quasar luminosity, at the very 
highest luminosities. Any quasars more luminous than My = — 27 should be found in 
massive elliptical hosts which at the present day would have My — —24.5. 

Key words: quasars: general - galaxies: active ~ galaxies: evolution - black hole 
physics 



1 INTRODUCTION 

Thanks largely to the resolution advantage offered by 
the Hubble Space Telescope (HST) the last decade 
has seen huge advances in our understanding of 
the host galaxi es of th e nearest (z < 0.3) q uasars 

jBahcall ct al] Il997l iHooper. Irnpev fc Foltd Il997l 

Bovcc ct al. 199^; McLurc ct al. 1999; McLcod fc McLeoc 
200L, .Harnilton. Casertano fc Turnshck 2002 , Dunlop et al. 
2003l) . HST observations have demonstrated that low-z 
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quasar hosts are luminous {L > L*) galaxies, confirming the 
results of earlier ground-based studies. But the key advan- 
tage of HST has been its ability to distinguish between disc 
and elliptical morphologies, leading to the finding that all 
radio-loud quasars (RLQs) and the majority of radio-quiet 
quasars (RQQs) reside in massive bulge-dominated galaxies. 

With recent improvements in the capabilities of both 
HST and ground-based telescopes, new studies are be- 
ginning to shed light on the evolution of quasar hosts 
from high redshifts (z > 1) down to the present 
day ||Lehner^^^J_^99|, j^lomo^^otilain«i^&TVeved 
I2OOII IStockton fc Ridgwavl I2001I: iRidgwav et alj I2OOII: 
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iKukula et al.l200ll:lH'utchings et al.l2002l) . At the same time 
it has become increasingly clear from studies of inactive 
galaxies that black-hole and galaxy formation and growth 
are intimately linked processes, resulting in the now well- 
established correlation between black-hol e mass and the 
mass of the host galaxy's stellar bulge l|M ^rorri m et alJ 
Il998l: iGebhardt et"ai]l2000l : iFerrarese fc Merrittil200(]r 



Most previous studies of quasar hosts have concentrated 
on quasars of relatively low luminosity, largely because it 
is much easier to disentangle the host and nuclear light in 
such objects. However, the quasar population spans lumi- 
nosities ranging from the (admittedly somewhat arbitrary) 
transition from Seyfert galaxies at My ~ —23 through 
to the most luminous objects, with absolute magnitudes 
Mv ~ —30; a factor of ~ 1000 in terms of luminosity. The 
majority of quasars currently known at large redshifts be- 
long to the bright end of the luminosity function. This is due 
to the degeneracy between redshift and luminosity in flux- 
limited samples, a situation which is beginning to be rec- 
tified as modern deep surveys detect low-luminosity, high- 
z quasars in increasing numbers. However, to understand 
the behaviour of the most massive galaxies and their black- 
holes in the cosmologically-interesting high redshift regime 
{z > 2) will inevitably require the study of the most lumi- 
nous quasars at these redshifts. 

The aim of the current study is to help break the degen- 
eracy between quasar luminosity and redshift by studying a 
sample of quasars at a single redshift that spans an apprecia- 
ble fraction of the quasar luminosity range (see Fig.l). The 
lowest redshift at which this can be attempted is z ~ 0.4, 
since the cosmological volume enclosed at lower redshift is 
too small to contain examples of the rare, most luminous 
quasars, with My < —27. Not only does the design of this 
study allow us to explore the relation between quasar lumi- 
nosity and host-galaxy properties, but the most luminous 
objects in this programme can also provide a low-redshift 
baseline against which to compare the hosts of luminous 
high-z quasars in future studies. 

The key technical difficulty in this work is, of course, 
the disentanglement of the point spread function (PSF) 
produced by the bright unresolved nucleus, from the emis- 
si on produced by th e host galaxy. In this study, as 
in lDunlop et al.l ll2003ll . we have used 2D modelling to per- 
form this decomposition. 

The paper is structured as follows. In Section 2 we de- 
scribe the quasar sample along with our selection criteria. 
Section 3 details our observing strategy and Section 4 de- 
scribes the reduction of the HST images. The 2-D modelling 
procedure used to analyse the images and extract informa- 
tion about the hosts is described in Section 5. Section 6 sum- 
marises the results produced by this image modelling, with 
the images themselves being located in an Appendix, along 
with detailed information about each object. In Section 7 
we discuss our principal results and their implications for 
our understanding of the quasar phenomenon. Finally, our 
conclusions are summarised in Section 8. 

For ease of comparison with our previous work we adopt 
an Einstein-de Sitter universe with Ho = 50 km s~^ Mpc~^. 
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Figure 1. Absolute magnitude versus redshift for quasars ob- 
served to date in our HST host-galaxy imaging programmes. 
Filled circles represent radio-quiet quasars (RQQs) and open cir- 
cles radio-loud quasars (RLQs). Our earlier work (small symbols) 
concentrated mainly on quasars of moderate luminosity (typi- 
cally Mv > -25) in three redshift regimes (z ~ 0.2, 1 & 2), 
allowing us to probe the evol ution of the host galaxies over a 
large fraction of cosmic history iMcLure et al .ll99alKukula et alj 
2001; DunloD et al. 2003). The current study (large symbols) is 
designed to explore an orthogonal direction in the My — z plane, 
by sampling a large range of quasar luminosities at a single red- 
shift, 2 ~ 0.4. This is the lowest redshift at which one can find 
very luminous quasars (those with My < —27), comparable to 
the most luminous quasars in the high-redshift universe. 



2 THE QUASAR SAMPLE 

The sample was selected from the quasar catalogue of 
IVeron-Cettv fc VeronI Jl993h and comprises two subsamples, 
both confined to the redshift range 0.29 < z < 0.43 (Ta- 
ble 1). The first, 'low-luminosity' subsample consists of five 
radio-loud and five radio-quiet quasars (RLQs & RQQs) 
with absolute magnitudes —24 > My > —25. All of the 
RLQs have 5 GHz radio luminosities > 10'^'* W Hz~^sr~^ 
and steep radio spectra to ensure that their radio lumi- 
nosities have not been significantly boosted by relativistic 
beaming. The RQQs have all been surveyed in the radio at 
sufficient depth to ensure that their 5 GHz luminosities are 
indeed < 10^* W Hz~^sr~^. The second, 'high-luminosity' 
subsample consists of all known quasars in this redshift 
range with absolute magnitudes My < —26, and includes 
2 quasars with My ~ —28. These two samples allow us 
to explore an orthogonal direction in the optical luminos- 
ity - redshift plane, in contrast to our previous HST stud- 
ies of quasar hosts jMcLure et al.lll999l:lKukula et alJl200ll : 
IDuuIod et alJl2003r which concentrated on quasars of com- 
parably moderate luminosity {My > —25), but spanning a 
wide range in redshift out to 2; ~ 2 (Fig.l). 
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Table 1. Quasars in the current study. J2000 co-ordinates were obtained from the Digitised Sky Survey plates maintained by the Space 
Telescope Science Institute. Redshifts and apparent V magnitudes are from the quasar catalogue of Vcron-Ccttv Sz Vcron |200oj). For 
consistency we use a B1950 lAU format to refer to the quasars in this paper; the name under which each object appears in the HST 
Archive is given in the final column and additional names are given in the description of each object in the Appendix. The low-luminosity 
subsample was observed in HST Cycle 7, using the WF2 chip and the F814W filter. The high-luminosity subsample was observed in Cycle 
9 using the WF3 chip, and the slightly narrower F791W filter. This latter observing run included one additional object, 1404—049, an 
inactive spiral galaxy at redshift 0.04 which had been misclassified as a quasar. This object has been omitted from the analysis presented 
in this paper. 



Name (B1950) Type z RA (J2000) Dec (J2000) V My Observing date HST archive name 



Low-luminosity subsample (Cycle 7; F814W) 



1237-040 


RQQ 


0.371 


12 


: 39 


39.0 


-04 


16 


38.5 


16.96 


-24.77 


Feb 12 1999 


1239-041 


1313-014 


RQQ 


0.406 


13 


: 16 


09.0 


-01 


54 


54.1 


17.54 


-24.39 


Feb 01 1999 


Q1313-0138 


1258-015 


RQQ 


0.410 


12 


: 58 


13.9 


-02 


00 


09.3 


17.53 


-24.42 


Feb 19 1999 


1258-015 


1357-024 


RQQ 


0.418 


14 


: 00 


06.4 


-02 


42 


22.6 


17.43 


-24.56 


Mar 01 2000 


1400-024 


1254+021 


RQQ 


0.421 


12 


: 57 


05.9 


+01 


49 


46.8 


17.14 


-24.86 


Feb 11 1999 


1257+015 


1150+497 


RLQ 


0.334 


11 


: 53 


22.3 


+49 


30 


21.4 


17.10 


-24.40 


Nov 25 2000 


LB2136 


1233-240 


RLQ 


0.355 


12 


: 35 


39.6 


-25 


11 


31.1 


17.18 


-24.46 


Aug 09 1997 


PKS1233-24 


0110+297 


RLQ 


0.363 


01 


: 13 


22.0 


+29 


58 


58.8 


17.00 


-24.68 


Feb 02 1999 


B2-0110+29 


0812+020 


RLQ 


0.402 


08 


: 15 


21.3 


+01 


55 


44.8 


17.10 


-24.80 


Feb 18 1999 


PKS0812+02 


1058+110 


RLQ 


0.423 


11 


: 00 


49.3 


+10 


47 


00.3 


17.10 


-24.90 


Apr 02 1999 


AO1058+11 


GRW+70D5824 


Star 




13 


: 38 


58.0 


+70 


16 


32.1 


12.77 




Feb 25 2000 


PSF-STAR 










High-luminosity subsample 


(Cycle 9; F791W) 






0624+691 


RQQ 


0.370 


06 


: 30 


08.6 


+69 


05 


40.0 


14.20 


-27.73 


May 25 2000 


HS0624+6907 


1821+643 


RQQ 


0.297 


18 


: 22 


02.8 


+64 


20 


05.3 


14.10 


-27.31 


Aug 17 2000 


E1821+643 


1252+020 


RQQ 


0.345 


12 


: 55 


22.5 


+01 


43 


46.2 


15.48 


-26.29 


May 16 2000 


EQS-B1252+020 


1001+291 


RQQ 


0.330 


10 


: 04 


06.1 


+28 


55 


19.2 


15.50 


-26.16 


May 21 2000 


TON0028 


1230+097 


RQQ 


0.415 


12 


: 33 


28.8 


+09 


31 


04.9 


16.15 


-26.05 


Jul 5 2000 


1230+097 


0031-707 


RLQ 


0.363 


00 


: 33 


57.2 


-70 


25 


24.5 


15.50 


-26.39 


Aug 7 2000 


MC4 


1208+322 


RLQ 


0.388 


12 


: 10 


39.8 


+31 


56 


26.0 


16.00 


-26.04 


May 7 2000 


B2-1208+32 


GRW+70D5824 


Star 




13 


: 38 


54.1 


+70 


16 


21.2 


12.77 




May 1 2000 


PSF-STAR 



3 OBSERVING STRATEGY 

All of our previous HST observations of quasar hosts were 
carefully designed to maximize the chances of successfully 
separating the starlight of the host from the PSF of the 
central quasar. We used the same observing strategy for the 
current observations and, since some of the quasars in our 
new sample are significantly more luminous than those in the 
earlier programmes, these precautions assume even greater 
importance. 

Observing dates for each object are listed in Table 1, 
along with the name under which the dataset is listed in the 
HST Archive. The observations were carried out in two dif- 
ferent HST observing cycles, although in practice the dates 
overlap. Observations of the low-luminosity subsample were 
carried out in Cycle 7 whilst the high-luminosity subsample 
was observed as part of Cycle 9. 

3.1 Choice of Alter 

As in our previous programmes we selected filters to corre- 
spond to V-band in the quasar's restframe. This ensures that 
our images sample the object's restframe spectrum long- 
wards of the 4000 A break, where the starlight from the host 
is relatively bright, whilst avoiding strong emission lines such 
as Ha and [Olll]A5007. Despite being directly associated 
with the quasar activity, ionised emission-line regions can 
extend over several kiloparsec. By excluding such emission 
from the images, we obtain a cleaner picture of the distri- 
bution of starlight in the hosts. 



For the low- luminosity subsample we used the F814W 
'broad /' filter which corresponds closely to the standard 
Cousins /-band. The high-luminosity subsample spans a 
slightly broader range of redshifts and in order to avoid 
contamination of the images by emission lines we used the 
slightly narrower F791W filter. 

3.2 Choice of detector 

Observations were made with the HST's Wide Field & Plan- 
etary Camera 2 (WFPC2). We opted to use the WF chips 
(800 X 800 pixels, with a scale of 100 mas pixeP^) since 
their relatively large pixels offer greater sensitivity to low 
surface brightness emission. Targets were centred on one of 
the three WF chips, the exact choice depending on which of 
the three had performed best over the period immediately 
prior to the Phase 2 proposal deadline. 

3.3 Exposure times 

High dynamic range is imperative in a study of this kind 
in which it is necessary to accurately characterise both the 
central core of the quasar as well as the faint outer wings 
of the PSF and the underlying host. The approach used is 
the same as that tried and tested in our previous works (e.g. 
McLure et al. 1999, Kukula et al. 2001, Dunlop et al. 2003). 
We took several exposures of increasing length, with expo- 
sure times carefully scaled so that no image would saturate 
beyond the radius out to which the PSF could be followed in 



4 D. J. E. Floyd et al. 



Table 2. V-band magnitudes and 5 GHz radio flux densi- 
ties for the quasars in our sample. At the time of writing the 
majority of the RQQs have only upper limits to their radio 
fluxes, but these are sufficient to place them safely below the 
accepted radio-loud/radio-quict boundary of L^qjj^ < 10~^* 
W Hz~ ^ sr~^. The final column give s the sour ce of the radio data 
(VCV = IVeron-Cettv fc VeronI <200(t: G-l-99= lGoldschmidt etld] 
il99Hh : B+96= lBlundell et al.l il99fl) . Upper limits from the 
NVSS were converted to 5 GHz by assuming a radio spectral 
index a = 0.5 (/^ oc v~^-^). 



Source 


Class 


z 


V 


SsGHz 
(mjy) 




Ref. 






Low- 


-luminosity 








1237-040 


RQQ 


0.371 


16.96 


< 0.3 


< 


22.14 


G+99 


1313-014 


RQQ 


0.406 


17.54 


< 0.3 


< 


22.21 


G+99 


1258-015 


RQQ 


0.410 


17.53 


< 0.2 


< 


22.05 


G+99 


1357-024 


RQQ 


0.418 


17.43 


< 0.3 


< 


22.24 


G+99 


1254+021 


RQQ 


0.421 


17.14 


< 0.3 


< 


22.25 


G+99 


1150+497 


RLQ 


0.334 


17.10 


717.0 




25.62 


VCV 


1233-240 


RLQ 


0.355 


17.18 


670.0 




25.45 


VCV 


0110+297 


RLQ 


0.363 


17.00 


340.0 




25.17 


VCV 


0812+020 


RLQ 


0.402 


17.10 


845.0 




25.62 


VCV 


1058+110 


RLQ 


0.420 


17.10 


225.0 




25.12 


VCV 






High 


-luminosity 








0624+691 


RQQ 


0.370 


14.20 


< 1.2 


< 


22.75 


NVSS 


1821+643 


RQQ 


0.297 


14.10 


8.6 




23.58 


B+96 


1252+020 


RQQ 


0.345 


15.48 


0.8 




22.50 


G+99 


1001+291 


RQQ 


0.330 


15.50 


< 1.2 


< 


22.65 


NVSS 


1230+097 


RQQ 


0.415 


16.15 


< 1.2 


< 


22.85 


NVSS 


0031-707 


RLQ 


0.363 


15.50 


95.0 




24.62 


VCV 


1208+322 


RLQ 


0.388 


16.00 


91.0 




24.90 


VCV 



the previous, shorter exposure. The series of exposures were 
then spliced together in annuli to construct an unsaturated 
image of the requisite depth (the pointing stability of HST 
between successive exposures using the FGS fine tracking 
mode is ~ 0.003 arcsec). 

For the 'low-luminosity' subsample a single orbit of HST 
time was sufficient for each object, with exposures of 5, 26 
and 3 x 600 seconds. For the more luminous quasars in 
the second subsample we required some shorter exposures 
to avoid saturation of the central pixels, as well as more 
long exposures to provide greater depth, since the wings of 
the quasar PSF encroach further out into the surrounding 
galaxy. Here we devoted two orbits to each object, with ex- 
posures of 2, 26, 2 X 100 and 3 x 600 seconds in the first 
orbit and 3 x 700 seconds in the second. 



3.4 PSF determination 

Our previous host-galaxy studies have emphasized the im- 
portance of characterising the instrumental PSF over a large 
range in angular radius in order to accurately separate the 
contributions of host galaxy and active nucleus. The struc- 
ture of the HST PSF is quite variable, especially at large 
radii, and depends on its position in the instrument field of 
view, the SED of the target source and the timing of the 
observations. 

We therefore devoted two orbits of our HST time 
allocation to obtaining deep, unsaturated stellar PSFs 
through both the F814W and F791W filters. The star used 



was GRW+70D5824 {V = 12.77), the same white dwarf 
used in our previous quasar host-galaxy studies with HST 
iMcLure et al .ll999l:lDunlop et alJl2003ll . This star is an op- 
tical standard for WFPC2, and so very accurate photometry 
is readily available. Its DA3 spectral type and neutral colour 
{B — V = 0.09) mimics well the typical quasar SED at the 
redshift of our sample. In addition, as there are no com- 
parably bright stars within 30 arcsec, we can be sure that 
our stellar PSF is not contaminated by light from nearby ob- 
jects. All observations were performed using the same region 
of the WF chip, within the central 50x50 pixel region. 

For the PSF observations we used an observing strategy 
similar to that for the quasars in order to obtain final images 
with high dynamic range. A series of exposures were carried 
out with durations of 0.23, 2.0, 26.0 and 160.0 seconds. Al- 
though the increasing exposure times lead to increasing sat- 
uration in the core, they never become saturated outside the 
radius at which the signal-to-noise of the preceding (shorter) 
exposure has become unsatisfactory. A 2-point dither pat- 
tern was also adopted in order to achieve better sampling 
(the /-band PSF is under-sampled by the 0.1-arcsec pixels 
of the WF chips). Such a strategy is not employed for the 
quasar observations, since the read noise introduced by the 
additional observations results in an unacceptable decrease 
in signal-to-noise ratio. 



4 DATA REDUCTION 

All the images were passed through the standard WFPC2 
pipeline software which performs many of the initial im- 
age processing and calibration steps such as dark and bias 
frame subtraction, along with flat fielding. We carried out 
three additional reduction steps prior to analysing the im- 
ages; cosmic-ray decontamination, sky subtraction, and re- 
construction of the saturated core regions of both quasar and 
stellar PSFs by splicing in images with shorter exposures. 



4.1 Removal of cosmic rays and bad pixels 

Cosmic rays were removed using CRREJ, an iterative 
sigma-clipping algorithm in iraf, which rejects high pix- 
els from sets of exposures of the same field. The result is a 
single deep image (with integration time equal to the sum 
of its parts). This technique however cannot remove the nu- 
merous persistent bad or "hot" pixels that appear in the 
WFPC2 chips. These are listed in a tabular form on the 
WFPC2 web page and can be incorporated into a mask for 
modelling purposes (see Section 5). 



4.2 Background subtraction 

In order to remove a smoothly varying sky background from 
each image we adopted the following procedure. First, the 
mean and standard deviation of the background in each im- 
age were estimated from a subset of pixels excluding obvious 
sources. An estimate of the sky distribution was then created 
by fitting a 2nd order polynomial to each image, having re- 
placed la deviations away from the mean background level 
by the mean value. This smooth sky map was subtracted 
from each original image to give a final sky-subtracted frame. 
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The background on WFPC2 images is extremely flat, and 
higher order approximations found to be unnecessary. 



4.3 Building deep, unsaturated quasar images 

After cosmic ray removal from the deep images, the central 
regions from the shorter exposures were cut out, multiplied 
by a factor determined by annular photometry, and sub- 
stituted into the saturated cores. All steps were performed 
using standard data-reduction tasks within iraf. For the 
brightest quasars in our sample, in which the saturated re- 
gion was large, we spliced together three images in this way 
to create the final image. 



5 IMAGE ANALYSIS AND MODELLING 

The techniques and software used to perform the 
2D modelling have been discus sed previously in detail 
jMcLure. Dunlop fc Kukulal l2000 h . and in this Section we 
only briefly reconsider the procedure, highlighting the mod- 
ifications required to cope with the extreme luminosities of 
some of the quasars in the current study. 



5.1 Quasar models 

Each model quasar was constructed by combining a model 
host galaxy with a central delta function to represent the 
nuclear point source. We model the host galaxy surface flux 
using the Sersic equation tSersic ,1968,l : 



E(r) = Eo exp 



.Ro) 



(1) 



The surface brightness E(r) describes an azimuthally- 
symmetric distribution, which is projected on to a gen- 
eralised elliptical coordinate system to allow for different 
eccentricities and orientations of the host galaxy. Rq is 
the characteristic scale length, Eo is the central surface 
brightness, and /3 describes the overall shape of the pro- 
file. We follow e d the procedure of iMcLure et al.l il999l) and 
iDunlop et al.l (|200^ in fitting first to a priori elliptical 
(/3 — 0.25) and disc (/3 = 1.0) models, and then to a more 
general form in which /? is a free parameter. Throughout 
this work we use the term "scale length" to refer to the 
half-light radius of the host galaxy models; i.e. the radius 
within which half of the integrated light is contained. For 
elliptical galaxy models, this half-light radius is equivalent 
to the often quoted effective radius (or Re), while for the 
disc or spiral models {j3 — 1.0) the half-light radius and the 
exponential radius are related by R1/2 = 1.687?o- 

Care must be taken in translation from the continuous 
distribution described by equation 1, to the discrete regime 
of CCD pixels. Simply basing the surface brightness of each 
pixel upon the radius of the centre of that pixel is inaccu- 
rate, particularly toward the centre, where E(r) varies non- 
linearly across the angular width of a pixel. To account for 
this, galaxy models were calculated at a much higher resolu- 
tion, and then rebinned to WF chip resolution for compari- 
son with the data. Each pixel value depends upon at least 9 
calculations of surface brightness, and up to 676 calculations 



within the central 0.5 arcsec. Flux is then added to the cen- 
tral pixel of the model galaxy, to represent the unresolved 
nuclear component. 

A given model quasar is thus specified uniquely by 
a point in the 6-dimensional parameter space {X} — 
{Ljv,Eo,ro,/3, 6,e}: 

• Ln = luminosity of the nucleus 

• Eo = central surface brightness of the host galaxy 

• ro = characteristic scale length of the host galaxy 

• = position angle of the host galaxy 

• e = J = axial ratio of the host galaxy 

• (3 controlling the shape of the profile 

The result is an idealised, seeing- and diffraction-free image 
of a quasar, which can then be convolved with the appropri- 
ate PSF to produce a simulated HST observation. 



5.2 Modelling the PSF 

HST/WFPC2 offers both extremely deep imaging and an ex- 
tremely well characterised PSF. However, the PSF is under- 
sampled on the WF chips and care was therefore taken to 
match the central regions of the PSF to each quasar image. 
The sub-pixel centring of each quasar image was found us- 
ing the CENTROID routine in iraf to characterise the distri- 
bution of light in the central region. Accurate oversampled 
models of the central regions of th e PSF were then gen- 
erated using the tinytim software jKristlll999ll and were 
re-sampled using the correct central position to provide an 
accurate model of the central few pixels of each quasar im- 
age. 

The TINYTIM calculation depends upon optical path dif- 
ferences within the Optical Telescope Assembly, and can 
be performed for any sampling rate, or position within the 
WFPC2 field-of-view. Using 21-times oversampled (with re- 
spect to a WF pixel) PSFs rebinned to the WF chip resolu- 
tion, we matched the sub-pixel position of the centre of the 
PSF to each quasar through a 2D minimum grid search. 
The best-fit tinytim model was then scaled up (by annular 
photometry) and spliced into the centre of the deep stel- 
lar PSF image. This PSF could then be convolved with the 
model galaxy plus nuclear component and the result com- 
pared to the real quasar image. 



5.3 Pixel Error Analysis 

The error allocation for such a technique must be done care- 
fully if the figure of merit is to have real meaning. Inac- 
curate error weighting may lead to the dominance of one 
region over another in the fitting, and hence to biased re- 
sults. Pixel values were assumed to be independent and to 
obey Poisson statistics. We used a combination of Poisson 
and sampling errors. 



5.3.1 Poisson Noise Error 

The minimum possible noise in a given pixel is the combina- 
tion of Poisson error due to photon shot noise, dark current, 
and the noise introduced by CCD read-out. However, this 
simple error calculation severely underestimates the effec- 
tive error in the central region of our quasar images, due to 
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sub-critical sampling of the steeply rising PSF by the WF 
pixels. 

5.3.2 Sampling Errors 

For each quasar image we constructed a "PSF residual" 
frame by scaling up the PSF to match the quasar in the 
central pixel and subtracting. The resulting frame can be 
used to quantify the extent to which the re-sampled tiny- 
TlM-|-stellar PSF constructed for that particular quasar im- 
age has succeeded in mimicking the central, pixelised flux 
distribution. Specifically, we compute the inferred remain- 
ing sampling error as a function of radius by calculating the 
variance of the distribution of pixel values in one-pixel wide 
annuli in this residual image. This variance (cr^) value is then 
assigned to each pixel in a given annulus. This is done from 
the centre out, until the value of this sampling error has 
fallen (as expected at some radius, given adequate knowl- 
edge of the PSF) to the level of the mean Poisson noise error 
discussed above. We call this radius the "sampling error ra- 
dius", Rsamp. We could choose to use this approach across 
the whole image, since it falls to the same level as the aver- 
age Poisson noise by a radius of typically ~ 1 arcsec from the 
centre of the quasar image. However, it is clearly preferable, 
where possible, to assign each pixel an error based upon its 
own individual noise properties, rather than a blanket error 
for an entire annulus. The contribution of genuine galaxy 
residuals to this error calculation are negligible. 

For the two most luminous objects in this study, E1821 
and HS0624, the sampling error in the PSF residual image 
remained above the mean Poisson noise at radii much larger 
than 1 arcsec. This is because, in these two cases the central 
quasar is so bright that the image actually contains more 
information on the detailed structure of the PSF at large 
radii than does our deepest image of the PSF star. Conse- 
quently, for these two objects, the errors in our knowledge 
of the PSF at radii of several arcsec become important, and 
we had to enhance the adopted errors in the image at large 
radii (by typically \/2; in practice we adopted an average of 
the 'sampling' and Poisson errors at large radius) to achieve 
an acceptable model fit with a flat distribution of values in 
the final image produced by the model-fitting procedure. 

5.3.3 Central pixel 

No sampling error can be computed for the single central 
pixel. In this case we apply the Poisson error, noting that the 
central pixel value has been scaled up from a very short (0.26 
seconds) snapshot exposure in order to avoid saturation. The 
error on the central pixel is typically of the same order as, 
or a little larger than, the sampling error deduced for the 
innermost annulus. 

5.4 Minimization 

The model quasar was convolved with the PSF and, after 
masking of diffraction spikes and nearby companion objects, 
compared with the HST data. The minimum was found 
for each object within the 6-dime nsional parameter space us- 
ing the downhill simplex method dPress et alJll993) . Six dif- 
ferent starting points (spanning the entire parameter space) 



were used for each object to ensure that the global mini- 
mum had been located. This model was checked using three 
starting simplexes at ±2, 5, and 10% of each best-fitting 
parameter value. Finally, the model was double-checked by 
performing a grid search around the best-fit model. 



5.5 Photometry 

Photometric calibration was performed using the HST head- 
ers PHOTFLAM and PHOTZPT, in order to convert 
counts into physical units of spectroscopic flux density (erg 
s^^ cm~^ ^~^)- In each case, the rest frame fllter band is 
calculated and compared to standard Johnson V-band. We 
assume a simple flat spectrum across the fllter bandpass. 
The internal uncertainty in photometry due to the accu- 
racy of the calibration reference files and the stability of the 
instrument is 1-2%. The conversion from the HST to the 
Johnson photometric system also has an uncertainty of a 
few percent. 



6 RESULTS 

As described in the previous section, we used two separate 
modelling strategies in order to determine the morphology 
of the hosts and the relative contributions of the nuclear 
and galaxy components. In the first case we fitted a pure 
de Vaucouleurs (r^'^^-law) elliptical galaxy and then a pure 
(exponential) Freeman disc to the host and used the differ- 
ence in the x^ values for the two models to decide which 
model gave the best overall fit. Unless otherwise stated in 
the notes, all objects were modelled out to a radius of 4 
arcsec. Table 3 lists the results of this strategy. 

In the second case, we carried out modelling using a 
variable- /3 fit, in which the /3 parameter of equation 1 is 
allowed to vary freely. This allows for a more general mor- 
phology than the strictly disc or bulge technique. Table 4 
shows the results of this variable-/? fitting which, for the 
most part, rather impressively reinforces the results of the 
fixed-/3 models. However there are a few objects in which the 
variable- /9 technique returned a hybrid value and these are 
noted in the entry for the relevant object in the Appendix. 

Greyscale images of the individual objects are also pre- 
sented in the Appendix. For each quasar we show the final 
reduced /-band (F814W/F791W) HST image (top left), the 
best-fitting model quasar (nuclear point source plus either 
pure bulge or pure disc host galaxy) to the quasar image 
(top right), the model host galaxy only (bottom left) and 
the model-subtracted residual image (bottom right). Radial 
profiles for the best-fit bulge and disc models are presented 
in Fig.2. 

We can gain further insight into how successful we 
have been in disentangling the host galaxy from the nu- 
cleus through investigation of the contours in the — 
Ri/2 plane (e.g. Fig. 3). For any quasar in which we have 
successfully characterised the host luminosity (i.e. elimi- 
nated the degeneracy between host and nuclear contribu- 
t ion s), these contours should lie along a s l ope of 5 .0 (se e 
e.g. lAbraham. Crawford fc McHardvl fl99l iMalkanI ITgsll . 
and allow us to assess how well constrained these two pa- 
rameters are. 
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Table 3. Results of model fitting with both de Vaucouleurs spheroid and Freeman disc models. Columns are as follows: object name; 
best fitting host-galaxy morphology (disc or elliptical); reduced-x^ value for the best fit model; Ax^ between the chosen and alternative- 
morphology model; half-light radius, -R1/21 of best fitting galaxy model in kpc; surface brightness of the host at the half-light radius, 
A'l/2i ™ units of V mag.arcsec"^; integrated absolute magnitudes of the nucleus and the host galaxy, converted from the appropriate 
filter band (F814W/F791W) to Johnson \/-band; the ratio of integrated nuclear and host galaxy luminosities; position angle of the host 
(in degrees, anti-clockwiso from vertical in the images); the axial ratio of the host. 



lAU name 


Morphology 
(best fit) 




Ax^ 


Rl/2 

(kpc) 


Ml/2 




J^,fhost 


Ln/Lh 


PA 

(°) 


a/b 










Radio-Quiet Quasars 












0624-f691 


Elliptical 


1.488 


486.2 


9.7 + 0.7 


22.0 + 0.15 


-27.18 


-24.01 


18.45 


140 


1.2 


1001+291 


Elliptical 


1.318 


326.3 


15.4 + 0.6 


23.1 + 0.10 


—25.62 


—23.47 


7.27 


58 


1.7 


1230+097 


Elliptical 


1.219 


107.7 


5.8 + 0.2 


21.7 + 0.13 


-25.24 


-23.24 


6.26 


1 


1.3 


1237-040 


Disc 


1.323 


334.5 


6.7 + 0.1 


22.0 + 0.05 


-23.98 


-22.63 


3.46 


14 


1.1 


1252+020 


Elliptical 


1.228 


102.0 


3.9 + 0.5 


22.1 + 0.30 


-25.26 


-22.04 


19.47 


150 


1.1 


1254+021 


Elliptical 


1.356 


646.3 


14.2 + 0.3 


23.2 + 0.05 


-23.91 


-24.00 


0.92 


30 


1.1 


1258-015 


Elliptical 


1.352 


14.4 


1.5 + 0.2 


19.8 + 0.25 


-23.77 


-22.33 


3.77 


140 


1.1 


1313-014 


Disc 


1.254 


351.0 


5.6 + 0.1 


21.7 + 0.06 


-23.74 


-22.68 


2.65 


174 


1.1 


1357-024 


Disc 


1.257 


457.1 


5.8 + 0.1 


21.9 + 0.05 


-23.66 


-22.47 


2.99 


160 


1.2 


1821+643 


Elliptical 


1.828 


570.4 


18.9 + 0.2 


22.9 + 0.05 


-27.14 


-24.33 


13.35 


113 


1.3 










Radio-Loud Quasars 












0031-707 


Elliptical 


1.268 


802.8 


11.0 + 0.4 


23.1 + 0.08 


-23.84 


-23.21 


1.80 


70 


1.2 


0110+297 


Elliptical 


1.327 


561.1 


12.3 + 0.8 


23.6 + 0.12 


-23.93 


-22.99 


2.39 


31 


1.2 


0812+020 


Elliptical 


1.509 


904.7 


17.4 + 0.3 


23.6 + 0.05 


-24.80 


-23.81 


2.49 


155 


1.2 


1058+110 


Elliptical 


1.390 


135.5 


13.1 ± 1.1 


24.1 + 0.14 


-23.58 


-22.69 


2.28 


159 


1.3 


1150+497 


Elliptical 


1.360 


340.6 


8.3 + 0.3 


22.1 + 0.11 


-24.09 


-23.28 


2.11 


174 


1.5 


1208+322 


Elliptical 


1.096 


49.1 


6.5 + 0.1 


22.2 + 0.05 


-25.01 


-22.50 


10.06 


7 


1.9 


1233-240 


Elliptical 


1.358 


47.7 


3.1 + 0.1 


20.6 + 0.05 


-24.78 


-22.98 


5.23 


58 


1.1 



Table 4. Outcome of variable-/? modelling. Columns are as fol- 
lows: object name; best-fit morphology from pure bulge & disc 
models (Table 3); best-fit /3 value with no assumed morphology; 
the value of reduced-x^ produced by this best-fit model; improve- 
ment in fit, Ax^ obtained by using the variable-/? technique com- 
pared to the best-fit fixed morphology model. 



lAU name 


Morph. 


/3 




Ax^ 




Radio-Quiet 1 


Quasars 






0624+691 


Elliptical 


0.20 


1.485 


42.2 


1001+291 


Elliptical 


0.26 


1.318 


0.1 


1230+097 


Elliptical 


0.37 


1.216 


12.0 


1237-040 


Disc 


0.96 


1.321 


0.6 


1252+020 


Elliptical 


0.22 


1.227 


0.7 


1254+021 


Elliptical 


0.24 


1.356 


2.0 


1258-015 


Elliptical 


0.26 


1.351 


2.0 


1313-014 


Disc 


0.97 


1.254 


0.2 


1357-024 


Disc 


1.32 


1.246 


22.8 


1821+643 


Elliptical 


0.22 


1.771 


485.4 




Radio-Loud ( 


Quasars 






0031-707 


Elliptical 


0.26 


1.266 


11.9 


0110+297 


Elliptical 


0.22 


1.323 


15.5 


0812+020 


Elliptical 


0.23 


1.503 


9.5 


1058+110 


Elliptical 


0.33 


1.388 


6.9 


1150+497 


Elliptical 


0.36 


1.356 


19.4 


1208+322 


Elliptical 


0.32 


1.091 


19.5 


1233-240 


Intermediate 


0.56 


1.361 


57.2 



7 DISCUSSION 

The quasars imaged in this study span almost two orders 
of magnitude in optical luminosity but only a narrow range 



of redshifts. They therefore allow us to investigate the rela- 
tionship between galaxies and their central black-holes, and 
the relative roles of black-hole mass and fuelling efficiency 
in determining quasar luminosities. 

We have successfully recovered a host galaxy for each 
one of the 17-strong sample. In general, the host size and 
central surface brightness have been constrained to within 
a few kiloparsecs, and half a magnitude, respectively, as is 
illustrated by the typical joint confidence region illustrated 
in the left-hand panel of Fig. 3. However, there are two ob- 
jects, the RQQ's 1252+020, and 1258-015 for which the fit 
yields poor stability for the host properties, and the result- 
ing much larger confidence region for 1252+020 is shown 
in the right-hand panel of Fig. 3. Overall host and nuclear 
fluxes are typically constrained to better than O.lmag by 
the modelling software, with a similar uncertainty due to 
the conversion from ST mags to standard V^-band. 

7.1 Host galaxy morphologies 

With regard to basic host galaxy morphology, the results 
of this study arc quite clear cut, and confirm and extend 
the findings of DunloD ct al. (2003). For every quasar host 
the modelling software yielded a clear decision in favour of 
either a disc-dominated or bulge-dominated host. Moreover, 
in virtually every case this preference was confirmed by the 
variable-/? model, which returned a value of (3 very close to 
either 0.25 (elliptical) or 1 (disc). 

At this point it is important to clarify what we mean by 
"bulge-dominated" or "disc-dominated" galaxies. Our mod- 
elling software finds the best overall fit to the light from the 
quasar and its host galaxy. This is dominated by contribu- 





r/arcaac r/arcsac r/accaec 




(g) 1258-015 (E.QQ): EUiptwial, (h) 1313-014 (EQQ}: Dis*;, (i) 1357-024 (KQQ): Disc, 

J?.j^2 = 1.5kpfi 5.61q)c; ^1/2= 5.flkp(; 

Figure 2. The radial profiles of the best-fitting bulge or disc models for the 17 quasars in our sample. Each plot shows the azimuthally 
averaged image data (open circles with lu error bars), the aaimuthally averaged best- fit model after convolution with the PSF (solid 
line) and the azimuthally averaged best-fit unresolved nuclear component after convolution with the PSF (dotted line) . The form of the 
fit (disc or elliptical) and the scale length {R1/2) of the model galaxy are also given beneath each panel. 



tions from the point-like nucleus itself, and from the smooth, 
high SNR host region far from the nucleus. Our HST images 
are of sufficient depth that we expect to bo able to detect 
large features at least as dim as = 27 mag.arcsec^^. To 
place this in some context, this is a sufficient depth that the 
prominent tidal arm in Mrkl014, could be easily be detected 
if it were placed at a ^ ~ 1. ^ However, close to the nucleus, 



this sensitivity is impaired by our lack of knowledge of the 
form of the PSF, and here a small-scale, relatively bright 
feature might go unnoticed by the model. Such features can 
be exposed in our residual images, which show the best fit 
model subtracted from the data. Thus when we claim to 
detect bulge or disc-dominated hosts, we mean just that: 
on large scales of a few kpc, the host light is dominated by 



^ Maxkaxian 1014 is a well-known disturbed active galaxy at at a surface brightness of ~ 24 V mag. arcsec ^ (see images and 

z = 0.163, with a prominent tidal arm that is easily detected profiles in McLure et al. 1999). 
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(j) 1821+643 (R.QQ): EUiptifial, (k) 0031-707 (KLQ): EUipdkal, 



J?, 



:i^2 ^ ll.Okpc 



(1) 0110+297 (KLQ): ElUptical, 
Rif2 = 12.31q)c; 




(m) 0S12+O20 (ELQ): EUiptkal, 
= 17.4kp(! 



(n) lata +110 (KLQ): Ellirriikal, 
J?.i^2 = 13.1kp(; 



(o) llfiO+497 (KLQ): ElUptkial, 
Ai^2 =S.3kpfi 




(p) 1208+322 (KLQ): Elliptical, (q) 1233-240 (KLQ): Elliptical, 



R. 



li^a -6.5kpc 



Figure 2. - continued 



smooth emission that follows a spheroidal (r^^^'-iaw) or a 
disc-like (exponential) profile. There obviously remains the 
possibility of a small, centrally concentrated bulge compo- 
nent in the disc-dominated hosts, and in an effort to address 
this, we attempted a 2-component modelling technique. This 
is simply an extension of the modelling technique described 
above, with individual disc and bulge components each de- 
scribed in full by 4 parameters, giving a 9-dimensional model 
overall. The results are presented in table 5: A significant 
spheroidal component is found in two of the three disc- 
dominated hosts; 1237-040 and 1313-014. 

As illustrated in Fig. 4, the three quasars in the cur- 
rent sample for which we find disc-dominated hosts are, i) 



radio-quiet, and ii) in the low-luminosity sub-sample with 
Mv > —25. In fact, reference to Table 3 reveals that the 
three disc-dominated hosts house nuclei with My > —24. 
Furthermore, the brightest of the three quasars (1237—040), 
is found to have a significant spheroid. This result therefore 
meshes well with the luminosity-depe ndence of host - galax y 
morphology illustrated in Fig. 10 of iDunlop et al.l i2003l) : 
disc-dominated host galaxies are not found for nuclei with 
Mv < —24. In one other case (1313—014) we detect a 
spheroidal component at low luminosity {My = —20.28), 
suggesting an accretion ra te around L/LEdd — 3, as suming 
the local MBH-Msph ratio lMcLure fc Dunlonl J2002l) . In the 
final disc-dominated object (1357—024) we detect no signif- 
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M 0,8 0,85 0,5 0,55 0,S 0,65 

logio Ri/gApc login ^i/g/kpc 



Figure 3. x' contours (at 1 — Sc levels) in the fi — R plane for 1150+497 (left), and 1252+020 (right). In each case, the best-fit solution 
is marked by a dot. The figure demonstrates the degeneracy that remains between galaxy size and surface brightness, even when any 
confusion between host and nuclear light has been eliminated. 1150+497 is a typical example of the objects in this study and the contours 
describe a slope close to 5 (dashed line), as expected if the host galaxy's luminosity has been correctly constrained. Of all the objects 
in our sample, 1252+020 has the least robust host luminosity constraint, the contours lying along a somewhat steeper slope of ~ 5.7. 
One other object, 1258—015, exhibits slightly steep contours, having a slope of ~ 5.6. In these latter 2 cases there is clearly still some 
degeneracy between host and nuclear flux. 



Table 5. Results of the 2-componcnt model fitting to the disc-dominated hosts. The improvement in the fit over the simple disc case is 
given in column 3. The nuclear luminosity, and the scale-length and luminosity of each host galaxy component are given in columns 4-8, 
with the bulge: disc ratio in the final column. 



Source 


Morph. 


Ax' 


Mv(Nuc) 




My (Bulge) 




Mv (Disc) 


Bulge/Disc 










kpc 




kpc 






1237-040 


D/B 


10.8 


-23.96 


3.8 


-21.18 


6.5 


-22.46 


0.31 


1313-014 


D/B 


4.7 


-23.72 


3.3 


-20.28 


5.6 


-22.62 


0.12 


1357-024 


D 


0.3 


-23.65 


0.4 


-18.38 


6.3 


-22.47 


0.023 



leant bulge eomponent. However, for a roughly Eddington- 
limited accretion rate, we would anticipate a bulge luminos- 
ity of around Mv = —20.8. Such a bulge is expected to have 
a scalelength of ~ 3kpc, and could be difficult to distinguish 
from the unresolved nuclear light. 

As discussed by iDunlop et al.l (|200^, this result can 
now be viewed as a natural consequence of the now well- 
established proportionality of black-hole and spheroid mass. 

7.2 Host galaxy scale lengths and luminosities 

Table 3 lists the scale lengths and luminosities for the best- 
fit fixed-morphology models. O nce again, the result s are 
broadly consi s tent w ith those of iMcLure et alJ lll999l) and 
iDunlop et al.l i2003ft : the hosts are generally large, luminous 
galax;ies. 

Three of the five smallest galaxies are the disc- 
dominated hosts. There is a tendency for the hosts of the 
RLQs to be slightly larger than those of the RQQs, but this 
is not statistically significant. 



{i?l/2)(flQQ) =8.7±1.8kpc 

On average the more luminous quasars also have slightly 
larger hosts, but again the mean values for the two subsam- 
ples are in agreement given the statistical uncertainty. 

{i?l/2)(A/„<-25) = 10.0±2.4kpc 

(-Ri/2)(A/jv>-25) = 9.0 ± 1.5kpc 

We find that that all the hosts are more lumino us than 
L* {M^ = -21.0: lEfstathiou. Ellis fc Petersonll988h . There 
is no statistically significant difference between the average 
values for each subsample, but these basic statistics should 
not obscure the fact that the two quasars in the sample with 
Mv{Nuc) < —27 are also the only two objects for which we 
find Mv{Host) < -24. 

{Mv{Hoat))^nLQ) = -23.06 ± 0.16 

{Mv{Host))(iiQQ) = -23.12 ± 0.25 

{Mv{Host))^Mr,<-25) = -23.27 ± 0.36 



{Rl/2)(RLQ) = 10.2±1.8kpc 



{Mv{Host))^M^>-25) = -23.01 ± 0.16 
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RQQ 



RLQ 



Mv>-25 



Mv<-25 



Figure 4. Histograms of the best-fit /3 values from our variable- 
P models, f) = 0.25 is equivalent to an r^/** de Vaucouleurs law, 
/3 = 1 is an exponential Freeman disc. In the panel to the left 
we divide the sample in terms of radio luminosity. On the right, 
the sample is divided according to optical absolute magnitude. 
Clearly all optically-luminous, and all radio-loud objects lie in 
bulge-dominated hosts, confirming and extending the trends de- 
duced by Dunlop ct al. (2003). 



7.3 Kormendy relation 

The Kormendy relation is the photometric projection of the 
fundamental plane exhibited by elliptical galaxies. The host 
galaxies of the quasars in our sample follow a Kormendy 
relation of the form 



^1/2 = (19.2 ± 0.6) + (3.33 ± 0.7)logioRi/2 



(2) 



shown in Fig. 5 (where we have plotted and fitted only 
those with bulge-dominated hosts). A galaxy with a well- 
constrained luminosity but unknown scale length will lie 
along a locus with a slope of 5, illustrated by the error ellipse 
in the top right corner of this figure (c.f. Fig. 3). The slope of 
3.33 is in excellent agreement with that determined recently 
for 9000 early-type galaxies (3 .33 ± 0.09) drawn from the 
SDSS bv lBernardi et aO i2003l) and is sufficiently different 
to a slope of 5 to convince us that the surface brightnesses 
and scale lengths of the hosts have been well constrained. 



T.4 The role of galaxy mergers and interactions 

Interactions and merging events between galaxies have long 
been suggested as the triggering events for the activation of 
quasars, especially at low redshifts where some mechanism 
is required to initiate fuelling of the black-holes in otherwise 
stable, gas-depleted ellipticals. Indeed, most host galaxy 




0.0 0.5 1.0 1.5 
logio Ri/a/kpC 

Figure 5. The scale length {R1/2) vs surface brightness (Ati/2) 
projection of the fundamental plane. Filled circles represent RQQ 
hosts, open circles RLQ hosts (bulge-dominated hosts only). The 
solid line shows the best fit Kormendy relation to the sample and 
has the form ^iya = (19-2 ± 0.6) -I- (3.33 ± Q.7)logioRi/2- The 
narrow ellipse in the top right corner of the plot shows the 2cr 
error contours for 1258—015; its slope of 5 is due to the remain- 
ing degeneracy between size and surface brightness when host 
luminosity has been well constrained. 



studies to date have found that indications of disturbance 
such as tidal tails, multiple nuclei and close co mpanions are 
present in around 50% of quasar ho sts (e.g . ISmith et alj 
[l986; Hutchings fc Ncff 1992; Bahcall et alJ Il997l. How- 
ever iDunlop et alJ ll2003l) point out that this is also true 
of inactive massive ellipticals, so that it is not clear whether 
mergers are genuinely a defining feature of quasar hosts 
or merely the legacy of their parent population. Certainly 
many quasar hosts appear to be entirely undisturbed, so 
clearly a large-scale disruption of the host is not always nec- 
essary to trigger fuelling of the central engine (or at least 
the timescales for relaxation after the merger event and of 
fuel reaching the central engine may sometimes be vastly 
different). 

The residual images of the quasars in the current study 
(see Appendix), provide a means of identifying signs of 
galaxy interactions which might not be obvious in the raw 
HST images. They are produced by subtracting the best- 
fitting axially-symmetric quasar model from the HST image. 
Since the model only attempts to fit the smooth underlying 
distribution of galaxy light, any additional structures (spiral 
arms, bars, tidal tails, double nuclei etc) will be made more 
obvious in the residual image. 

In our sample we find unambiguous evidence for an 
ongoing galaxy interaction in only one object, the RQQ 
1237—040. Several other objects are candidates for some 
form of disturbance having taken place (for instance the 
RQQ 1001-1-291 with its prominent spiral arms and large- 
scale de Vaucouleurs profile), or have other objects within a 
few arcsec on the sky which might conceivably be interact- 
ing if they lie at the same redshift. In fact the majority of 
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our quasars appear to have some companions nearby on the 
sky, which is at least suggestive of a cluster environment. 

However we find no correlation between the luminosity 
of the quasar and the presence of any morphological distur- 
bance in the host. In our small sample, at least, the most 
luminous quasars seem no more likely to be interacting sys- 
tems than their less luminous counterparts. 

7.5 Black- hole masses 

Reliable black-h ole masses are available for at least 37 
nearby galaxies jKormendv fc Gebhardtl |2001^ ■ highlight- 
ing the correlations which exi st between spheroid lum inos- 
ity and black-hole mass (e.g. 'Mago rrian et al.]ll998^ . and 
between spheroid velocity-dispersion and black-h ole mass 
jGebhardt et al.ll2000l: [Perrarese fc Merritti l200(tl . As dis- 
cussed b'v TMcLur^^Dunlopril2o'o2l) . both these relations for 
inactive ellipticals, and the results of H/3-derived virial black- 
hole estimates in active objects, are consistent with a direct 
proportionality of black-hole and spheroid mass of the form 
Mbh ~ 0.0012Afspfi , with a typical scatter of 0.3 dex. 

We used the luminosities from our best-fit bulge- 
dominated galaxy models to estimate host-galaxy masses, 
given an estimate of the mass-to light ratio of a n 
early- type galaxy jjgrgensen. Franx fc KiaergaardI Il99d) : 
{M/L)ii-band oc L^'^^ . We then used the black-hole:spheroid 
mass relation above to estimate black-hole masses for the 
quasars, independent of their observed nuclear output. The 
results of this calculation are given in column 3 of Table 
6. Once again, this has only been performed for the bulge- 
dominated hosts. A similar calculation can be performed for 
the spheroidal components of the disc-dominated hosts (see 
end of section 7.1). The median black-hole masses for the 
high and low-luminosity subsamples are: 

{Mbh){M„<-25) = (1-1 ± 1.0) X 10^ 

{Mbh)(m^>~25} = (0.5 ±0.3) X 10^ 

We find that all of the host galaxies are sufficiently mas- 
sive {Msph > IO^'^Mq) to contain a black-hole in excess of 
10* Mq, but the difference in mass between the black-holes 
in optically powerful and optically weak quasars is not large 
enough to account for the factor ~ 10 increase in luminos- 
ity. Unfortunately, no direct measures of black-hole mass are 
available for any of the objects in our sample. A compari- 
son with, for example, a Virially estimated black-hole mass 
would provide an extremely valuable test. 

7.6 Fuelling efficiencies 

We can now calculate the predicted luminosity of each bulge- 
dominated object if the black-hole were to radiate at its 
Eddington limit {L^f^ = 1.26 x 10^^ Watts) and com- 
pare this with the actual luminosity of the quasar nucleus 
obtained from our model fitting. The results of this proce- 
dure are listed in columns 4 and 5 of Table 6, and plotted in 
Fig. 6; there is clearly no correlation between black-hole mass 
and fuelling efficiency. Ifowever, one can immediately see a 
clearer distinction between our luminous and dim subsam- 
ples than was apparent simply from their estimated black- 
hole masses: 



Object Morph. Msph Mbh M^'^'' 

(lOiiAfe) (109Mq) 



Radio-Quiet Quasars 



0624-1-691 


E 


13.90 


1.67 


-27.28 


0.91 


1001-1-291 


E 


727 


0.87 


-26.57 


0.42 


1230-1-097 


E 


5.51 


0.66 


-26.27 


0.39 


1252-1-020 


E 


1.29 


0.16 


-24.70 


1.67 


1254-1-021 


E 


13.70 


1.65 


-27.27 


0.05 


1258-015 


E 


1.84 


0.22 


-25.08 


0.30 


1821-1-643 


E 


20.50 


2.46 


-27.70 


0.60 



Radio-Loud Quasars 



0031-707 


E 


5.31 


0.64 


-26.23 


0.11 


0110-1-297 


E 


4.07 


0.48 


-25.94 


0.16 


0812-1-020 


E 


11.00 


1.31 


-27.02 


0.13 


1058-1-110 


E 


2.84 


0.34 


-25.55 


0.16 


1150-1-497 


E 


5.78 


0.69 


-26.32 


0.13 


1208-1-322 


E 


2.26 


0.27 


-25.30 


0.76 


1233-240 


E 


4.02 


0.48 


-25.93 


0.35 



Table 6. Galaxy spheroid and black-hole mass estimates for each 
of the bulge-dominated quasars in our sample. The table also 
lists the theoretical Eddington luminosity, My'^'^, of each black- 
hole, and the efficiency at which the black-hole is accreting ex- 
pressed as the ratio of the luminosity ascribed by our model to 
the nuclear point source to the Eddington luminosity predicted 
by the model of the host galaxy (Lj^/LEdd)- Note that the RQQ 
1252-1-020 appears to be accreting at a super-Eddington rate. This 
object has the least robust model fit of the entire sample, and it 
is likely that the nuclear flux has been overestimated. Results for 
the spheroidal components of the three disc-dominated objects 
are presented where detected. Both detections (1237—040 and 
1313—014) yield slightly super-Eddington accretion rates, whilst 
the nuclear luminosity of 1357—024 can be explained in terms of 
Eddington-limited accretion in an undetected, compact central 
bulge at Mv ~ -20.8. 

(7^>(A/„<-25) =0.78 ±0.2 

^Edd 

<7^)(M«>-25) =0.17 ±0.03 

^Edd 

If we exclude the relatively poorly constrained luminous 
RQQ 1252±020 (which appears, from our modelling, to ex- 
ceed the Eddington limit), the median Eddington ratio for 
the luminous subsample drops to 0.47, but this is clearly 
still significantly higher than for the low-luminosity sample. 

Thus, within our z = 0.4 sample, increasing quasar lu- 
minosity appears, on average, to refiect a mix of both larger 
black-hole mass, and increased fuelling efficiency. Not sur- 
prisingly, the only two quasars in the present sample with 
Mv < —27 have the two most massive black-holes. 

A number of other features of the results summarised 
in Table 6 are worthy of comment. First, while inferred fu- 
elling efficiencies range over an order of magnitude, we find 
no evidence for super-Eddington accretion. If one excludes 
the poorly constrained RQQ 1252±020, the most efficient 
emitter is 0624±691, with Ljv/I/Bdd — 1. Second, the most 
massive central black-hole found in our sample has a mass 
of 3 X 10^ Mq, comparable to the inferred ma ss of the super- 
massive black-holes at the centres of M87 l*M arconi et alJ 
[l997) and Cygnus A ijTadhunter et al.l2 003). Thus, the basic 
physical quantities derived for the quasars in our sample ap- 
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Figure 6. Quasar efficiency as a function of the Eddington lumi- 
nosity versus black-liole mass, as determined from the spheroid 
luminosity. Only bulge-dominated objects are plotted, and they 
have been divided into optically luminous (open squares) and 
optically dim (solid squares) subsamples at My(Nuc) = —25, 
as discussed in the text. Our model fit to the RQQ 1252-f020 
implies that the object has a small host with an extremely lu- 
minous nucleus, yielding a super-Eddington luminosity. However, 
the quality of the model fit is poor, and this value is quite likely 
erroneous. Overall there is no obvious tendency for fuelling effi- 
ciency to vary as a function of black-hole mass, and the increase 
in luminosity across the sample reflects an increase both in black- 
hole mass, and accretion efficiency. 



pear to be entirely reasonable, requiring neither unorthodox 
methods of accretion, nor surprisingly massive black-holes. 

It is important to note that the black-hole mass cal- 
culation applied above (and hence the values of Mbh and 
Ln/LEdd given in table 6 and Fig. 6) assume a single, fixed 
value for the black-hole:spheroid mass ratio. At some level, 
this is clearly unrealistic and it is therefore not in fact obvi- 
ous to what extent the scatter in the nuclear:host ratio re- 
flects a range of fuelling efficiency. Accordingly, we conclude 
this paper with the first exploration of the extent to which 
scatter in nuclear:host ratio can or cannot be explained by 
intrinsic scatter in the underlying black-hole:spheroid mass 
relationship. 



7.7 Black- hole mass versus fuelling rate 

The nature of the link between quasar luminosity and black- 
hole mass is more easily explored by plotting host versus 
nuclear luminosity. This is shown in Fig. 7 where we have 
plotted the absolute magnitudes of the hosts against those 
of the nuclei in our sample (circles), with 100%, 10% and 1% 
of the Eddington limit shown as solid, dash ed and dotted 
lines r espectively. Show n also are points fro mlDunlop et al 
('2003') (diamonds) and iMcLeod. Rieke fc Storrie-Lombard: 



Q999 ) (triangles), converted to rest-frame l^-band, and our 
adopted c osmology. We h a ve als o included 3 objects from the 
sample of iPercival et al.) ||200 j) (stars), for which archival 
HST images are now available (0043-1-039, 0316-346 and 
1216+069). It now seems likely that seeing limitations in 



this ground-based study effectively prevented successful dis- 
entanglement of host and nuclear fluxes, and accurate mor- 
phological distinctions. The replacement images from the 
HST archive have been analysed in precisely the same way 
as has been described for the present sample, and converted 
into rest-frame l/-band. 

The top panel of Fig. 7 shows that, while central black- 
holes appear to accrete with a wide range of efficiencies, the 
objects we term quasars are generally produced by black- 
holes emitting at > 10% of their Eddington limit, residing 
in host galaxies with L > L* . However, perhaps the most im- 
pressive feature of this plot is that, for a given host galaxy lu- 
minosity, the most luminous nuclear source has a luminosity 
essentially exactly as would be predicted from the Eddington 
limit corresponding to the mass of the central black-hole as 
deduced from the relationship Mbh ~ 0.0012Msph. In other 
words, while the statistical correlation between host-galaxy 
and nuclear luminosity within these samples may not be very 
strong, the relationship between host-galaxy and maximum 
nuclear luminosity appears extremely tight, and completely 
consistent with Eddington lim ited accretion. Su c h a re sult 
has previously been noted by iMcLeod fc Riek3 1^9^ for 
lower-luminosity AGN, but to the best of our knowledge, 
this is the first time that it has been demonstrated concretely 
at high luminosity. So clean is this relation over two orders 
of magnitude in Ljv that it has the potential to constrain the 
size of the scatter in the underlying black-hole:spheroid rela- 
tion for massive galaxies. In turn, such constraints can then 
illuminate the extent to which the apparent 1 dex scatter in 
fuelling efficiency can also be explained by intrinsic scatter 
in the underlying black-hole:spheroid mass relationship. 

A full exploration of this approach is beyond the scope 
of the present paper and is in any case better applied to 
larger and more statistically complete samples. However, we 
have produced the lower panel of Fig. 7 to illustrate the ex- 
tent to which the data can be reproduced by folding in an 
assumed underlying scatter in the black-hole;spheroid mass 
ratio, while assuming a single value for the Eddington ratio. 

We generated a random population of spheroidal galax- 
ies, and central black-holes defined by a Schechter function 
with a = 1.25 and /offio( *^|°' ) — 11.5 (where Mtot denotes 
the turnover in the distribution of total mass, stellar and 
dark matter), and assuming a fixed scatter in the black- 
hole:spheroid mass relation. From this galaxy / black-hole 
population we generate a sample of quasars, assuming a 
fixed fuelling efficiency, and constrained to share the same 
nuclear luminosity distribution as is found for the real sam- 
ple in the top panel of Fig. 7. 

Adoption of a scatter larger than 0.3 dex produces sig- 
nificantly more apparently super-Eddington objects than 
are observed. Conversely, adoption of a scatter substantially 
smaller than 0.3 dex can reproduce the apparently tight Ed- 
dington limit more closely, but underpredicts the apparent 
scatter in fuelling efficiency, and cannot therefore reproduce 
the tight Eddington limit displayed by the data in the top 
panel of Fig. 7. 

We used a 2D K-S test to compare the simulation with 
the data. By marginalising over universal efficiency we find 
the probability, p, that the quasar sample is consistent with 
a fixed fuelling efficiency. We find that for a 0.3dex or smaller 
scatter, the quasar sample is inconsistent with a population 
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Figure 7. Upper panel: Host versus nuclear luminosity for 
the bulge-dominated quasars in the current sample (circles). 
The spheroidal components of the two disc-dominated quasars, 
1258—015 and 1237—040 are shown by asterisks. The samples of 
McLeod et al. (1999) (triangles), and Dunlop et al. (2003) (dia- 
monds), plus three objects from Percival et al. (2001) re-imaged 
with the HST (stars - see main text) are also presented. Filled 
symbols once again denote radio-quiet objects, and open symbols 
radio-loud objects. The solid, dashed and dotted lines represent 
objects radiating at 100%, 10% and 1% of the Eddington lumi- 
nosity respectively on the assumption of a fixed black-hole to 
bulge mass ratio of Mbh = 0.0012Msph. The majority of the 
quasars in our sample appear to be radiating at >10% of their 
Eddington limit (the single super-Eddington object is the poorly 
constrained 1252+020 - see Fig. 3.). Error bars are smaller than 
the symbols for all objects except for 1252-f020 and 1258-015. 
Lower panel: The predicted distribution on the host versus nu- 
clear luminosity plane for a sample of quasars radiating at 50% 
of the Eddington limit, given an adopted scatter of 0.3 dex in 
the black-hole:spheroid mass relation, and including the effect of 
the exponential cutoff in the luminosity function above L* (dot- 
dashed line). The random sample has been resampled in order 
to reflect the the same distribution of nuclear luminosities as is 
displayed by the combined data in the upper panel. This scenario 
re-produces much (but not all of) the observed scatter in apparent 
fuelling efficiency without in fact requiring a range of Eddington 
ratios. However, at the same time an assumed scatter of 0.3 dex or 



in which there is a fixed fuelling efficiency (p = 0.022). In- 
creasing the scatter produces an improvement to acceptable 
levels at 0.4dex (p = 0.157) and 0.45dex(p = 0.198), but a 
scatter as large as O.Sdex is strongly excluded (p — 0.002). 

The figure presented illustrates the situation for a sam- 
ple of quasars radiating at 50% of the Eddington limit, with 
an assumed scatter in the underlying black-hole:spheroid 
mass relation of 0.3 dex. These values were chosen for this 
illustration as the combination which best reproduces both 
the apparently tight Eddington limit, and level of scatter dis- 
played by the data in the upper panel of Fig. 7, adequately 
describing the spread in nuclear luminosities at My < — 23. 
However, for dimmer quasars, we find that a variable fu- 
elling rate is essential in order to explain the full range in 
■r^ observed. 

In summary, from simulations of the sort described 
above, the observed apparent tight upper (Eddington) limit 
on fuelling efficiency can be used to set an upper limit of 
0.3 dex on the scatter in the underlying black-hole:spheroid 
mass relation, consistent with other recently der ived values 
iMcLure fc Dunlodl2002l: iMarconi fc Huntll2003l) . A signif- 
icant fraction of the scatter observed in the upper panel of 
Fig. 7 can then still be attributed to the scatter in the un- 
derlying mass relationship, but some variation in assumed 
efficiency would still seem to be required to reproduce the 
most underluminous objects. 

Finally, we can use the tightness of the bound placed by 
the Eddington line to place constraints on the hosts of higher 
redshift quasars, where direct imaging of the host is not 
possible. A quick inspection of Fig. 7 reveals that any quasar 
brighter than Mv = —27 must be found in, or at least end 
up within, a spheroidal galaxy brighter than My = —24, and 
must be radiating at a rate close to its Eddington limit. The 
turnover in the Schechter function places a natural limit on 
the abundance of such large galaxies, and the data appear to 
show just such a cutoff at a host luminosity of AIv — —24.5. 



8 SUMMARY 

Through the careful analysis of deep HST images, we have 
succeeded in determining the basic properties of the host 
galaxies of quasars spanning a factor of ~ 20 in luminosity, 
but within a narrow redshift range at z ~ 0.4. The sam- 
ple under study contains both radio-loud and radio-quiet 
quasars, and includes some of the most luminous quasars 
known in the low-redshift universe. 

Our results confirm and extend the trends uncovered 
in ou r previous HST-based studies of lower lu minosity ob- 
jects llMcLure et al.lll999l lOunlop et alJl2003h . Specifically 
we find that the hosts of all the radio-loud quasars, and 
all the radio-quiet quasars with My < —24 are giant ellip- 
tical galaxies, with luminosities L > L* , and scale lengths 
Re ~ 10 kpc. Moreover, the Kormendy relation displayed by 
these host galaxies is indistinguishable from that displayed 
by nearby, inactive ellipticals. 

From the luminosities of their hosts we have estimated 
the masses of the black-holes which power the quasars using 
the relationship Mbh = 0.0012Mspft, and hence, via compar- 
ison with the quasar nuclear luminosities, also the efficiency 
with which each black-hole is emitting relative to the Ed- 
dington limit. We find that the order-of-magnitude increase 
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in nuclear luminosity across our sample can be explained 
by an increase in characteristic black-hole mass by a factor 
~ 3, coupled with a comparable increase in typical black- 
hole fuelling efficiency. However, we find no evidence for 
super-Eddington accretion, and the largest inferred black- 
hole mass in our sample is Mbh ~ 3 x 1O^M0, comparable 
to the mass of the black-holes at the centres of M87 and 
Cygnus A. 

We explore whether intrinsic scatter in the underly- 
ing Mbh '■ Msph relation (rather than a wide range in 
fuelling efficiency) can explain the observed scatter in the 
Mv{host) : Mv{nuc) plane occupied by quasars. We find 
that the observed tight upper limit on the relation between 
Mv {host) and maximum Mv{'nuc), (consistent with the Ed- 
dington limit inferred from a single constant of proportion- 
ality in the Mbh ■ Msph) constrains the scatter in the 
underlying black-hole:spheroid mass relation to be 0.3 dex 
or smaller, but that this mass-relation scatter can indeed 
explain a substantial fraction of the apparent range in fu- 
elling efficiency displayed by the quasars, particularly those 
at Mv < —23. The scatter also explains objects with ex- 
ceptionally high nuclear-to-host luminosity ratios without 
the need for super-Eddington accretion rates. Finally, our 
results imply that, due to the cutoff in the Schechter func- 
tion, any quasar more luminous than My — —27 must be 
destined to end up in a present-day massive elliptical with 
Mv ~ -24.5. 
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APPENDIX A: IMAGES & NOTES ON 
INDIVIDUAL OBJECTS 

For each quasar we show the final reduced 7-band 
(F814W/F791W) HST image (top left), the best-fit model 
(pure bulge or pure disc plus nuclear component) to the 
quasar image (top right), the model host galaxy after re- 
moval of the nucleus (bottom left) and the model-subtracted 
residuals (bottom right). North and East are marked on the 
object frame (with the arrow pointing North). The object 
and model frames are contoured at the same surface bright- 
ness levels. The level of the lowest contour is given (in V 
mag.arcsec"^) at the top right of the object frame. Succes- 
sive contour levels are separated by 1 mag.arcsec"^ . 

Al Radio-Quiet Quasars 

0624-1-691 (HS0624-I-6907). One of the brightest quasars in 
the sky, this object has b een the subject of a c omprehensive 
multi- wavelength study jReimers et al.l^l995^ ■ which clas- 
sified the host galaxy as a massive elliptical, fitting a de 
Vaucouleurs profile with a nominal scale length of 1.8kpc to 
their PSF-subtracted image. The quasar appears to lie in a 
cluster, with a number of small companion objects visible in 
the field. 

We find the host to be best fit by a giant elliptical galaxy 
(7?i/2 = 10 kpc), with an extremely strong nuclear compo- 
nent {Ln/Lh = 18). The variable-/? fit returns a value of 
/3 = 0.20, again consistent with a pure de Vaucouleurs ellip- 
tical host. 

Due to the extreme luminosity of this quasar the 
masking for the diffraction spikes was applied over a much 
larger area than for the majority of objects in this study. 
The host galaxy contribution is obvious out to a radius of 
at least 5 arcsec, and we used a fitting radius of 6 arcsec to 



ensure that all detectable host light was used to constrain 
the model parameters. 

1001+291 (TON0028, PC 1001-^292, 2MASSi 
J100402 5-I-285535). This object was studied in some 
detail bv lBovce. Disnev fc BleakenI lll999ll . who claimed two 
galactic nuclei; one 1.92 arcsec (14.6 kpc) to the south-west 
of the quasar nucleus, and the other 2.30 arcsec (15.9 
kpc) to the north-east. However, this claim appears to be 
a res ult of over-subtractio n of the nuclear point source, 
since iMArauez et al.l i200ll) showed that the host possesses 
prominent spiral arms and a bar which crosses the nucleus 
from north east to south west, although they were unable 
to fit a surface-brightness profile. 

In the current data we also find spiral arms and a 
nuclear bar which is clearly visible in the residual image of 
this object. However, we find that the surface light profile 
of the underlying smooth component is very well fitted by 
a de Vaucouleurs law {R1/2 — 15 kpc, /? = 0.26), suggesting 
a bulge-dominated host. 

1230-^097 (LBQS 1230+0947). We find the host galaxy of 
this quasar to be an elliptical with a scale length of R1/2 = 
6 kpc. There are a number of other objects in the same 
field with possible evidence for a tidal interaction with the 
nearest object to the north. 

Allowing for a variable value of /3 yields a slightly better 
fit with (3 — 0.37, suggesting a somewhat intermediate 
morphology. 

1237-040 (EQS B1237-0359). This object appears to be 
interacting with a companion galaxy to the north, and the 
residual image shows that a tidal tail has been induced in 
the quasar host itself. The host is best fitted by a large 
{Ri/2 = 4 kpc) disc galaxy, with the variable-/? modelling 
returns a value of l3 — 0.96. The residual image shows some 
excess nuclear fiux which has not been accounted for by 
the pure disc fit, and the 2-component modelling reveals 
a luminous spheroidal component (Bulge/Disc=0.31) of 
moderate size {R1/2 ~ 3.8kpc). The luminosity of the 
dominant disc is reduced by 0.2mag. This is the most 
luminous of the disc-hosted quasars studied here, and it is 
interesting to note that it is also the one with the most 
significant detection of a spheroidal host galaxy component. 

1252+020 (E QS B1252+020, HE 125 2+0200). This is 
a radio-quie t l|Goldschmidt et all Il999l) . X-ray detected 
I Voges_eLalJ [19991) qua sar with a strong UV excess 
I Goldschmidt et al.|[l993) . Of all the objects in the current 
sample, this quasar proved to be the hardest for which to 
achieve an unambiguous model fit to the host galaxy. In- 
deed, as is shown in Fig. 3, we were unable to constrain the 
luminosity of the host to the same extent as for the other 
quasars: the contours in the fi — R plane have a slope 
slightly steeper than 5. One nearby companion had to be 
masked out, along with the diffraction spikes, before mod- 
elling could be carried out. In addition, there is a faint region 
of nebulosity directly to the north of the quasar. 

Our best fit model has an elliptical host with 
Ri/2 ~ 4: kpc and the highest nuclear-to-host ratio in the 
sample, Ln/Lh = 19.5, although as has been stated, the 
host and nuclear flux have not been completely disentan- 
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gled, and there is a large error associated. The /3 parameter 
modelhng also favours an elliptical host, with /3 — 0.22. 

1254+021 (EQS B12 54+0206). This radio-quiet 
ijGoldschmidt et al.l Il99!j) quasar shows very smooth 
extended emission, with only weak diffraction spikes, and 
an absence of nearby companions. We find that the object 
is best fitted by a large (-R1/2 ~ 14 kpc) elliptical galaxy 
and a weak nuclear component. The variable-/? fit confirms 
the host morphology, returning a value of /3 = 0.24. No im- 
provement is made with the addition of a discy component. 

1258-015 (EQS B1255-0143, 2MASSi J1258152-015918). 
A highly compact object, with an almost stellar appearance 
in the IfST J-band image. However, there is sufficient galaxy 
light to model, and we find that the best-fit host is a small 
elliptical galaxy with a half-light radius of just 1.5 kpc. This 
is confirmed by the variable-/? model which returns a value 
of /3 = 0.26, with no appreciable improvement to the fit. 
Little excess flux remains in the residual image suggesting 
that the model has accurately accounted for all the host 
galaxy light, and the radial profile is a good fit to the data 
throughout. However, the contours in R1/2 — are 
slightly too steep to be convinced that we have completely 
resolved host from nucleus. 

1313-014 (Q1313-0138, EQS B1313-0138, LBQS 
1313—0138). This compact object was only modelled out 
to 3 arcsec, beyond which the SNR falls below 1. The 
nuclear component of this quasar is relatively weak, with 
no prominent diffraction spikes visible in the image. Hence, 
despite the small angular size of the host, the model fit 
to the galaxy is robust. A spiral feature is visible in the 
residual image, with possible evidence for a bar passing 
through the nucleus. We find the host to be best fitted 
by a disc model with R1/2 — 3 kpc, and this is supported 
by the variable-/? model which returns a best-fit value of 
/3 = 0.97. A low-level spheroidal component is detected by 
the 2-component modelling, implying a roughly Eddington 
accretion rate, although the improvement in the fit is 
marginal. 

1357-024 (EQS B1357-0227, 2MASSi J1400066-024131). 
There are several fainter objects in the field, suggesting that 
1357—024 might lie in a relatively rich cluster environment. 
Diffraction spikes from a nearby bright star are visible in 
the southeast quadrant of the image. However, the quasar 
itself is a compact source with no discernible diffraction 
spikes. All companion objects, and the majority of the 
southeast region of the image were masked out of the fit. In 
addition, we only modelled out to a radius of 3 arcsec, where 
we run into the background. The modelling software shows 
a strong preference for a disc-dominated host, although the 
variable-/? model returns the unusual value of /? = 1.32. 
Most likely this is due to the nearby stellar diffraction 
spike leading to a flatter profile (and hence higher /?). 
The residual image shows a small amount of excess flux 
in the nucleus. No improvement is obtained with the 
2-component modelling. However, for an Eddington limited 
nuclear source we would anticipate a bulge component at 
Mv — —20.8, which could easily be concealed in a compact 



(< Skpc) central bulge component. 

1821+643 (E1821+643, IRAS 18216+6418, 8C 1821+643). 
The brightest quasar in the current sample, this object has 
been extensively studied at many wavelengths. Extremely 
luminous in the infrared and also with a strong X-ray com- 
ponent, this was one of the first radio-quiet quasars to be 
studied in detail at ra dio wav elengths and is known to con- 
tain a small radio jet jBlundell fc R awlings 2001) . 

Although this quasar appears to lie in a rich field, 
most of the surrounding objects are believed to be part 
of a background cluster at a redshift z ~ 0.6. A previous 
study resolved the host galaxy, finding it to be large, fea- 
tureless and red, but fa iled to determine its morphology 
jHutchines fc Nefjll99ll ). In addition, the nucleus itself is 
unusually red, indicating the presence of large quantities 
of dust, tho ugh no discrete dust lanes have been observed. 
iMcLeod &: McLeodI 1120011) separated the host and nucleus 
in their Jf-band NICMOS imaging study, finding a luminous 
elliptical galaxy of magnitude Mr ~ —26.7, with a nuclear 
component with Mr = —29.2, (when converted into our 
cosmology) . 

Because of the prominent diffraction spikes a larger than 
usual region of the image was masked prior to modelling. 
However, since extended flux is clearly visible in the image 
out to a radius of at least 6 arcsec, we therefore used this 
as our fitting radius. The quasar is best modelled as a large 
elliptical host {R1/2 — 10 kpc), with a strong nuclear com- 
ponent (Lm / Lh ~ 11). The variable-/? model is good accord 
with this decision {/3 — 0.22), with a significant improvement 
in the fit. 

The residuals accentuate the nebulous artifact some 
4 arcsec east-southeast of the nucleus, and also appear to 
show a spiral-like feature wrapping around to the northeast. 
Due to the PSF sampling problems encountered in this 
highly luminous object, it remains unclear as to whether 
the latter is a genuine feature of the host galaxy, or simply 
a PSF artifact. 



A2 Radio-Loud Quasars 

0031-707 (MC4, 2MASSi J00 34052-702552). A radio- 
loud quasar llGreeorv et al.l ll994') originally identified as a 
Magellanic object due to its proximity to the galactic plane. 
There are a number of companion objects, suggesting that 
the object lies in fairly rich cluster environment, with the 
potential for interactions with nearby objects. The HST 
image shows a relatively weak nucleus (the model fit gives 
a nuclear/host ratio of Lm/Lh — 2). The host is best 
fit by an elliptical galaxy model with R1/2 ~ llkpc (the 
variable-/? model gives a value of /? = 0.26, with only a 
slight improvement in the quality of the fit). 

0110+297 (B2-0110+29, 4C 29.02, 2MASSi 
J0113242+295815). Malkaij Jl984l) attempted to re- 
solve the host galaxy of this quasar from the ground but 
was prevented from doing so by poor seeing. This quasar 
appears fairly compact in our HST image, with prominent 
diffraction spikes and a number of other objects nearby 
on the sky, including a well-resolved spiral galaxy some 4 
arcsec to the east. We found the best fitting host to be a 
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large elliptical galaxy with R1/2 — 12 kpc, confirmed by the 
variable-beta fit, which returns a best-fit value of /3 = 0.22. 
There is a small symmetrical circumnuclear artifact present 
in the residual image. 

0812-^020 (PKS0812-H02, 4C -F02.23) An earl y study of 
this object llWvckoff. Gehren fc Wehingeil [ToSlh measured 
the extent of the nebulosity surrounding the quasar and 
found it to have a dia meter of some 88 kpc. S ubsequent 
work was carried out bv lHutchings fc Nefi lll99Cf) . who fit- 
ted an elliptical galaxy profile to the host and obtained a 
scale length of around 12 kpc (converted to our cosmology). 
They also claimed to find evidence for a tidal interaction. 

This quasar lies in a crowded region of sky, and 
consequently a great deal of masking was required before 
modelling could be carried out. However, the host itself is 
relatively bright, and the preference is for a large elliptical 
galaxy with a scale length of ~ 17 kpc. The variable-/? 
modelling confirms the morphology of the host, returning 
a best-fit value of /3 = 0.23. We find some residual nuclear 
fiux, but no strong evidence for any disturbance or interac- 
tion in the host. A major improvement in the fit is obtained 
by adding a moderate disc component (-R1/2 =3.7kpc, 
Bulge/Disc=8.2). There is a 0.2mag drop in the luminosity 
of the spheroidal component, but the nuclear component is 
unchanged. 

1058-hllO (AO1058-m, PKS 1058+llC, 4C 10.30). There 
are a number of apparent compani on objects, a n d the 
cl ustering ampUtude was studie d bv lYee fc GreenI ( ^87 ) 
fc lGreen fc Ye3 Jl984 . However, iBlock fc StocktonI (Il99lli 
found thes e objects t o be at a different redshift from 
the quasar. iHutchin'g^ I^S^) detected extended nebulosity 
around this object, but was unable to fit a radial profile. 

Although the active nucleus appears to be relatively 
weak in our image, the host and nuclear components proved 
quite difficult to separate. However the model did converge 
on a large elliptical host, with i?i/2 = 13 kpc. The variable- 
/3 fit returns a slightly intermediate value of jS = 0.33. 
No signs of major disturbance are visible, although some 
circumnuclear flux remains in the residual image. Addition 
of a low-level disc (Bulge/Disc=16.0) results in a significant 
improvement in the fit, with no change in the properties of 
the dominant bulge model. 



viously been conjectur ed to be associated with the quasar 
jHutchines et al.lll988ll . Our image also shows two objects, 
roughly 2 arcsec to the north fc north west of the quasar 
which were masked out along with the diffraction spikes 
prior to modelling. 

Our modelling procedure shows a strong preference for 
an elliptical host galaxy, with the best-fit model having 
a scale length of R1/2 = 8 kpc. However, the variable-/? 
modelling returns a best-fit value of /? = 0.36, intermediate 
between pure bulge and disc morphologies. The reason for 
this discrepancy may be apparent in the residual image of 
the object which shows several regions of excess flux to the 
south of the nucleus. 



1208-^322 (B2-1208+32, 7C 1208+3213). This quasar 
was detected as a sof t X-ray source by Einstein 
( Puchnarcwicz et alJll992h . Optically, it appears to be a 
compact object with a strong nuclear component. We find 
an underlying elliptical host with a scale length i?i/2 = 6.5 
kpc. A slight improvement to the quality of the fit is ob- 
tained by allowing 13 to vary freely, giving a best-fit value of 
/3 = 0.32. 

Although there are no obvious signs of interaction, 
there do appear to be a number of small, faint com- 
panion objects surrounding the quasar. This is the only 
radio-loud object in the current sample whose accretion 
efficiency appears to come close to the Eddington limit 
{LN/LEdd = 0.76). 

1233-240 (PKS 1233-24. [HB89] 1232-249). 
IWvckoff et al.l il98lll found an extended nebulosity with 
a diameter of some 166 kp c surrounding this quasar. Th e 
object was also imag ed bv IVeron-Cettv fc WoltieJ Jl99Cl) . 
who found an elliptical host with magnitude My ~ —22.7, 
but were not able to provide a scale length. 

Despite the prominent diffraction spikes of the strong 
nuclear component, some galaxy light is clearly visible in 
our image, and there are also several other objects in the 
field. The best-fit host is an elliptical with a scale length of 
about 3kpc. Examination of the radial profile shows some 
excess fiux compared to the pure elliptical model and the 
variable-/? model returns a value of /? = 0.56 suggesting 
that a significant disc component is also present. 



1150-1-497 (LB2136, 4C 49.22). Several previous attempts 
have been made to detect the host galaxy of this Opti- 
cally Violent Variable (OVV) quasar. Malkan (1984), us- 
ing the Palomar 1.5m telescope, and seeing-degraded mod- 
els of elliptical and disc- like hosts, claimed to find a mas- 
sive elliptical host of scale length 25 kpc (when converted 
to our adopted cosmology). However, an exponential disc 
was f ound to give a reasonable fit by Hutchings ( 1987) fc 
iHutchings et alj lll98d). after PSF-subtraction, and ID pro- 
file fitting. Finallv lWright. McHardv fc Abrahainl lll99d) . by 
assuming an elliptical galaxy model (the relatively poor sam- 
pling in their data meant that no real morphological clas- 
sification could be performed), detected a host in JC-band, 
with Mk = -27.3 ± 0.6. 

The object appears to be elongated along a north- 
south axis in the current HST image. There are several 
fainter objects some 10 arcsec to the NE which have pre- 
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